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Abstract 

Supercooled hquids display fascinating properties upon cooling such as the emergence of dy- 
namic length scales. Different models strongly vary with respect to the choice of the elementary 
subsystems (CRR) as well as their mutual coupling. Here we show via computer simulations of 
a glass former that both ingredients can be identified via analysis of finite-size effects within the 
continuous-time random walk framework. The CRR already contain complete information about 
thermodynamics and diffusivity whereas the coupling determines structural relaxation and the 
emergence of dynamic length scales. 

PACS numbers: 61.20.Lc, 61.20.Ja, 64.70.pm 
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Understanding the complex dynamics of supercooled liquids down to the glass transition 
as reflected, e.g., by the emergence of dynamic length scales [l, 2| is still a highly controversal 



problem [3|, Qualitatively speaking the system decomposes into^oupled elementary units 
(sometimes denoted cooperatively rearranging regions (CRRs) [5]). Presently discussed 
models differ dramatically with respect to the nature of this decomposition. For example in 
the mosaic approach [6|-[9| the complexity is fully embedded in the properties of the CRRs, 
denoted mosaics, whereas in the other extreme of the kinetically constrained models (KCMs) 
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13| the complexity is exclusively is related to the coupling while the elementary units 



are just trivial two-state systems. 

Here we show that it is possible via computer simulations of a standard model glass 
former to obtain a clear-cut decomposition and thus to learn about the nature of the CRRs 
as well as the coupling between them. The key idea is to directly extract coupling effects 
from studying the size-dependence under periodic boundary conditions. The analysis is 
performed in the continuous-time random walk (CTRW) framework which for the flrst time 
is quantitatively extended to larger systems. Two conclusions can be drawn. Firstly, our 
results explain the fact that the structural relaxation time displays major flnite size effects 
whereas there are hardly any for the diffusion constant. Secondly, we present a minimum 
model of the glass transition which reflects the observed coupling effects. Albeit similar in 
spirit to the KCMs many important differences are present. 

We have simulated the binary mixture Lennard- Jones system (BMLJ) which can be 
regarded as a prototype glass-former [l4|. Simulations have been performed in the NVT 
ensemble. For the smaller system we have have used a slightly shorter cutoff of rc = 1.8 
15| . In previous work it has been shown that a system as small as A'mm = 65 particles 



with periodic boundary conditions basically displays very similar thermodynamic as well 
as diffusive properties as a macroscopic system in the ran ge o f temperatures, accessible to 
computer simulations (in particular between Tc and 2Tc) 16|. Interestingly, for the same 



temperature range = 40 is too small p/7[. Thus, a BMLJ system with ^ Nminjs close 



to a system reflecting the effects of just a single CRR (more precisely: two CRRs jlSl]). 



Due to the microscopic resolution computer simulations are very well suited to extract 
relevant information about the coupling effects. One route which has been chosen by Biroli 



et al. is to immobilize a large system except 



or a sphere of variable radius r and to study the 



dynamics in this sphere in dependence of r 19| . One faces the problem to extract from this 
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highly non-equihbrium situation the relevant equilibrium properties. In this work we analyze 
the size-dependence of the BML J system. From the perspective of a CRR the coupling effects 
are switched on when increasing the system size. Interestingly, the structural relaxation time 
21 1 as well as the length scale of dynamic heterogeneities display significant finite-size 
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effects 
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23|. 



A key observable to characterize the relaxation is the incoherent scattering function S{k, t) 
which is defined as 

= (exp(2fc[Af(t)]) (1) 

where Ar* denotes the single-particle displacement vector during time t. The contribution of 
structural relaxation to the decay of S{k,t) can be determined by referring to the inherent 
structure (IS) trajectory, obtained after minimizing every configuration j2J], refiecting the 
minima of the potential energy landscape (PEL) [25]. The resulting incoherent scattering 
function is denoted Sis{k,t). Qualitatively, this procedure corresponds to a removal of the 
vibrational contributions. Conceptually similar is the use of the metabasin (MB) trajectory 



261]. A MB consists of an appropriately chosen set of nearby IS. The k-dependent 



relaxation time, which refiects the time a particle on average needs to move the distance 

2n/k, is defined via r(A;) = dt Sis{k,t). 

In Fig{T]we show the wave vector dependent relaxation time T{k) based on Sjs{k, t) where 
;he large k limit is significantly beyond the maximum of the structure factor. As shown in 
271 ] the metabasin (MB) dynamics of a system with = Nmin can be described as a CTRW 

in configuration space (with waiting times tmb) where spatial and temporal properties are 



decoupled 



271 ]. As a consequence the k-dependence of r(/c) can be written as 
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where 
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D (xl/ {tmb), oc D{tIib). (3) 

Due to the presence of dynamic heterogeneities, the distribution of waiting times ^{tmb) is 
very broad so that Tq, ^ (tmb)- 

How to characterize the dynamics for large systems (A^ ^ A!mm)? In this limit the MB 
trajectory would violate one of the crucial assumptions of the CTRW approach because due 
to the dynamic heterogeneities successive MB transitions are often performed by identical 
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FIG. 1: The k-dependent relaxation time T(k) vs k^'^ for different system sizes at T=0.5. Included 
are linear fits which work very well for k < 8. From the two adjustable parameters it is possible to 
derive the values of the diffusion constant D and the structural relaxation time Tq,. 



particles [27l]. In particular, Eqj3]is no longer valid for large N. There exists, however, an 
elegant solution to this problem by defining local waiting times Tiocai on the single-particle 
level in real space 29I, 3o|- Specifically, one considers a trajectory which is time-averaged for 
a time somewhat longer than the vibrational time scale. Whenever the particle has moved 
a fixed distance d (here: d = 1/3 in LJ-units where the precise choice is irrelevant) one 
identifies a transition process. 

Validity of the CTRW-behavior in configuration space implies that the corresponding 
single-particle trajectories can be also described as a CTRW where Eql3]has to be modified by 
replacing tmb with Tiocai- One might expect that this CTRW- type single-particle dynamics 
should prevail also when increasing the system size. Indeed, as can be seen from FigJT]Eqj2] 
is applicable for all in the complete fc-range under consideration. In Figj2]the resulting 
values of D{N) and Ta{N) are shown for different temperatures. Note that the diffusion 
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FIG. 2: The size-dependence of the diffusion constant D and the structural relaxation time in 
relation to its value for N = 260. Whereas the diffusion constant only displays very minor finite 
size effects, dramatic effects are observed for the structural relaxation time. The inset contains a 
linear plot of the finite-size effects for Tq together an estimation of the coupling constant q. 

constant indeed only displays small finite size effects whereas the relaxation time varies by 
nearly one order of magnitude. At first glance this may seem somewhat counterintuitive 
because upon decreasing the system size major finite-size effects are observed for a small- 
scale (large k) observable, i.e. r„, but not for a large-scale (small k) observable, i.e. D. 

Finally, we analyse the A^-dependence of ^pijiocai) as shown in Fig. El The data are 
discussed in three different directions. First, for = Nmin and T = 0.5 we have also 
included ^{jmb) showing a very good agreement with ip^riocai), thereby supporting the 
mapping of the waiting times in configuration space to those in real space. The presence of 



a few extremely long waiting times {jiocai] 
only covers half of the system size Nmir 



reflects the fact that on average a MB transition 



18| . Thus, there is a finite probability that some 



particles remain immobile even longer than the longest MB waiting time. Second, in the 
inset of Figj3] we verify the relation in Eqj3] oc D{t1^^^i) within statistical uncertainty. 
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FIG. 3: The waiting time distribution v'('^«ocaO for T = 0.5 and T = 0.6 for different system sizes 
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as well as ^p{tmb) for T = 0.5 and N = 65. In the inset the CTRW relation Tq, oc D{t^^^^i) is 



checked by comparing different system sizes between = 65 and N = 1040 for T = 0.5. 

Third, one clearly sees that long residences of slow particles are strongly suppressed when 
approaching the large-size limit. This observation is particularly pronounced for the lower 
temperature where up to O(IO^) particles finite-size effects are present. One might be 
tempted to conclude that the additional fluctuations around a given subsystem increase the 
mobility and shifts all time scales to shorter times. This, however, disagrees with the very 
small finite-size effects of the diffusion constant which implies that the first moment of the 
waiting time distribution is basically unchanged. 

For a system with = Nmin it has been shown that the state of the system can be, to 



a reasonab 
transition 



e approximation, characterized by a rate F, denoting the probability for a MB 
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26| . To discuss the nature of the observed finite-size effects we assume for 



reasons of simplicity that this approximation is strictly valid. Furthermore we introduce 
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p{T) as the equilibrium distribution of rates. A finite width of p{T) is equivalent to the 
presence of dynamic heterogeneities. For = Nmin it turns out that F is strongly related 
to the MB energy j26|. Thus, the observed invariance of the thermodynamic properties, i.e. 
the trivial scaling of p{e) (e: MB energies) upon increasing the system size, strongly suggests 
that p{T) for the subsystem does not change either. This also follows from the previously 
reported observation that in terms of MB transitions the total system behaves as if it were a 
superposition of independent subsystems of size Nmin [Slj- As a consequence the remaining 
effect of the coupling is to give rise to fluctuations of the rate of the tagged subsystem due 
to relaxation processes in adjacent subsystems without modifying the overall distribution 
p{T) {passive exchange process). 

This coupling has one important consequence. During a very immobile period, charac- 
terized by a small value of F, it is very likely that the system acquires a new (and on average 
larger) value of F. This immediately explains the numerical observation that the long-time 
contribution to (p^Tiocai) disappears when the system size increases and thus the fluctuations 
become faster. 

Furthermore, these observations allow us to understand the nature of the flnite-size ef- 
fects of Sjs{k,t). In particular one needs to understand why flnite-size effects occur for r„ 
rather than D. Both observables are fundamentally different because (at least for spatially 
uncorrelated jumps) the diffusivity is fully determined by the properties of a single jump 
whereas for complete structural relaxation basically every particle has to move, requiring 
a large number of elementary relaxation processes. Thus, the presence of passive exchange 
processes is relevant for Tq, but not for D. A more formal description of the resulting dif- 
ferences can be achieved by using the relations D oc (F) and oc (F~^) where the average 
is over p{T) [s^. For reasons of simplicity we consider the case where the fluctuations of 
F are extremely fast. Naturally, these fluctuations leave the average rate (F) and thus D 
unchanged. For (F~^), however, F needs to be substituted by (F). Then one obtains 

Ta{N) = ((F)-i) = 1/(F) < (F-i) = T^{N = N^,n) (4) 

where the inequality is a strict mathematical statement if dynamic heterogeneities are 
present. Of course, the less pronounced the fluctuations the smaller the ratio Ta{N) / Ta{N = 

-^min ) • 

A simple visualization of the physical scenario is presented in Fig. |H (1) The local 
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FIG. 4: Schematic sketch of the physical scenario in supercooled liquids. After a MB transition 
in the central subsystem the rates in the adjacent subsystems may change. In a minimum model 
with probably q an adjacent subsystem selects a randomly new rate according to the Boltzmann 



probability and with probability 1 — q there is no effect 



261]. 



system of size A^mm fully contains the information about the thermodynamics as well as the 
diffusion constant (first moment of waiting time distribution). This system can be very well 
characterized in the MB framework (given p{e) and r(e)). Since, at least in the analyzed 
temperature range, A^mm is roughly temperature-independent no dynamic length scales are 
involved to characterize, e.g., the temperature dependence of the diffusion constant. (2) The 
macroscopic system can be regarded as a superposition of systems of size A^mm- The coupling 
is too weak to change the thermodynamics (basically, this was the criterion to define Nmin) 
but gives rise to a dynamic coupling between the subsystems. For the minimum model, 
sketched in Fig. H] one can directly estimate the coupling constant q in the limit of small q 
and not too large N/Nmin [ssl via the relation 



ra{N) 



N 



1 



(5) 

with /3jv/ ~ 0.4 and tJ^N = Nmin)/ (t)) ~ 30 for the present system as determined in Ref. 
[2?! in the CTRW context. Comparison with the numerical data yields q ~ 0.008, see Figj2j 
This clearly shows that efficient coupling does not occur on the time scale of MB transitions 
but for the present case on the time scale of 4rQ,(A^ = Nmin)- 

The resulting correlations between the mobilities of adjacent systems will lead to the 
emergence of temperature-dependent length scales of dynamic heterogeneities. These corre- 
lations do not show up via a spatial correlation of instantaneous rates but rather from the 
spatial correlation of the relaxation behavior during sufficiently long time intervals (e.g. of 
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the order of Tq), as exactly captured by the four-point correlation function 3J]. For a simple 
model system one can show |33| how the coupling modifies the value of and how dynamic 
correlations and thus dynamic length scales emerge. 

The aspect of coupling is rephrased in the KCMs in a schematic way. Thus it is not 
surprising that several similarities to atomic glass-formers exist (e.g. in the FA model the 
diffusion constant does not depend on the coupling and the elementary length scale, i.e. a 



spin, is temperature-independent 



12j). However, some important differences remain. (1) 



Whereas a system of size A^mm, representing approx. two elementary units, already reflects 
many features of the macroscopic system, the same number of spins in the KCM naturally 
does not display any relevant behavior. (2) In principle a subsystem can always relax without 
constraints by its neighbors (see also Ref.[13]) since the coupling is via passive exchange 
processes rather than yes/no-decisions. (3) Strong and fragile systems can be understood 
by the same coupling rules by just changing the properties of the potential energy landscape 
of the subsystems 26|| whereas for the KCMs different rules have to be formulated. (4) 
The experimentally observed correlation between thermodynamic and dynamic properties 
is kept whereas the thermodynamics is considered irrelevant in the KCM approach. (5) The 
coupling strength can be directly extracted from simulations of a atomistic glass-forming 
system, allowing a semi-quantitative mapping between system and model. 

In summary, the model, sketched above, can be, on the one hand, regarded as a general- 
ization of the KCM whereas, on the other hand, it is a minimum model, capturing a realistic 
version of the facilitation effect. In particular it is compatible with the thermodynamic and 
dynamic observations for the model glass former reported in this work. We just note that 
there is no obvious way to map the present results on the mosaic approach since the basic 
length scale, i.e. the scale of the mosaics, is strongly temperature dependent. 

We greatly acknowledge the financial support by the DFG (SFB 458). 
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FIG. 5: (Fig. SI) A model system with two subsystems; see text for details. 

Supporting Online Material 



Discussion of a simple but illuminating model 

To substantiate the findings, related to the couphng effect, we analyze a minimum model 
which leaves the thermodynamics as well as the diffusivity unchanged upon coupling. It 
consists of two subsystems with a coupling constant q, characterizing the strength of the 
coupling. The model is sketched in Fig. SI. The elementary subsystem contains two states. 
The slow escape rate from the lower state is F^, the faster escape rate from the upper state is 
Tf. The relevance of dynamical heterogeneities is characterized by Fj/F^. The degeneracy of 
the high-energy state is proportional to iff, that of the lower state to ips (v?s + V'/ = !)• The 
Boltzmann probabilities pi are proportional to (pi/Ti. For reasons of simplicity we choose 
Pf = Ps = 1/2 which implies (fi = Fi/(Fs + Tf). For this subsystem one directly obtains 
(r) = 2/(F, + F/) and r„ = {t^)/{2{t)) = (1/2)(1/F, + 1/F/). Now we consider a larger 
system, containing two of these subsystems which are dynamically coupled. A relaxation of, 
e.g., subsystem 1 has two consequences. First, subsystem 1 acquires a randomly new chosen 
state with probability ipi. This condition just reflects the independence of two subsequent 
states, populated by subsystem 1. Second, with probability q subsystem 2 randomly chooses 
a state with probability pi whereas with probability 1 — g no change occurs. The value of 
q determines the strength of the coupling (g = 0: no coupling, q = 1 maximum coupling). 
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FIG. 6: (Fig.S2)The dependence of Tq,/ (t) on T f/Tg. Included are theoretical predictions Tq,/ (r) = 
Tf/{ATs) and Ta/{T) = Tf/{16Ts) for the limit of large dynamic heterogeneities. 

This procedure guarantees that the Boltzmann distribution in one subsystem and the average 
waiting time does not change as a consequence of the couphng to dynamical processes in 
the other subsystem. 

In what follows we show via straightforward Monte Carlo simulations in combination 
with analytical expressions that (i) indeed decreases with increasing system size and that 
(ii) the coupling gives rise to dynamic length scales. The key observable, recorded during 
the Monte Carlo simulations are the different moments (r") of the waiting time distribution 
of a single subsystem, is defined as / dtSo{t) = (r^)/(2(r)) where So(t) denotes the 
probability that a subsystem does not perform a relaxation process during a time interval 
of length t. 

In Fig.S2 the dependence of r^/ (r) is shown as a function of Fj/Fs for different values of 
q. We checked that in agreement with the setup of the model the value of (r) oc 1/D does not 
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depend on q. Several interesting observations can be made. (1) Tq decreases with increasing 
coupling. This naturally reflects the general system size dependencies of Tq,, reported in the 
main text. (2) For large dynamic heterogeneities the q-dependence disappears as long as 
g > 0, indicating some degree of universality for large dynamic heterogeneities. (3) In this 
limit one obtains Ta/ {t) = Tf/{4:Ts) for g = and Ta/{T) = Tf/{16Ts) for q > 0. These 
relations can be easily rationalized. For q — and F/ ^ the relevant long-time decay 
of So{t) is given by (1/2) exp(-Fsf), implying ^ 1/(2F,). Together with (r) ^ 2/F/ the 
observed relation immediately follows. For q > a slow relaxation process is only observed 
if ioT t — both subsystems are in the slow state. If we define ni2{t) as the probability 
that none of the systems has performed a relaxation process until time t one naturally has 
nuit) = (1/4) exp(— 2F5t). The prefactor expresses the fact that only in one out of four 
cases both subsystems are initially slow. After the relaxation of, e.g., subsystem 1 and as a 
consequence oi (ff ':$> (fs subsystem 1 will typically become fast afterwards. The large number 
of subsequent relaxation processes will for 5 > finally give rise to an exchange of the slow 
to the fast state in subsystem 2. Thus, for qVf the time difference between the initial 

relaxation process in subsystem 1 and the first relaxation process in subsystem 2 is much 
smaller than l/F^. As a consequence one has So{t) ~ ni2(t), yielding Ta/ {t) = Tf/{16Tg) 
as indeed observed. 

At a particular time the rates in both subsystems are uncorrelated to each other. However, 
for g > both subsystems are dynamically correlated if observed for a finite time interval; see 
Fig. S3. In agreement with typical data for glass-forming systems one observes a maximum 
correlation for t ^ Ta. In analogy to above one expects for t > Ta that So{t) = ni2{t) = 
(1/4) exp(— 2rst) which agrees with the numerical data. 

Finite-size effects for small coupling constant q 

In generalization to above we consider two subsystems with a general distribution of 
rates piT). For one tagged system we write for the relaxation function (in discrete notation) 
'S'(^) — ^iPi9i{t) where pi denotes the probability that the system is in the initial state 
i and gi{t) the probability that the system in this state has not relaxed until time t. For 
an uncoupled system and the presence of simple rate processes one has gi{t) :— g'^{t) — 
exp(— Fji). Since the thermodynamics does not change upon coupling the values of the Pi 
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FIG. 7: (Fig.S3)In analogy to the four-point correlation functions the dynamic correlations for the 
present system can be characterized via ni2(t) — [5*0 (t)]^ where n\2{i) denotes the probability that 
the total system does not relax during a time interval of length t. The solid lines correspond to 
the theoretical expectation for long times. 

remain the same. Furthermore, we define Gi = J dtgi{t) and = YliPi^i- "^^^ index 
always denotes the uncoupled case. Note that 1/ (r) = (F) := J^Pi^i- 

For general q one has to keep track of pairs of adjacent states (see above). However, a 
dramatic simplification is possible in the limit of small q. Restricting oneself to the lowest 
order effect of the coupling one can systematically neglect the forward-backward reactions 
which are at least of order g^. Thus, one can employ a mean- field type approach where 
a tagged system experiences fluctuations which are uncorrelated to the resulting passive 
exchange processes. These fluctuations induce a reorganization of the tagged system with 
probability q per relaxation process (more specifically: MB transition). As a consequence 
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the rate of change is Q := q/ {t). Formally, this can be written as 26 1 



{d/dt)g,{t) = -Vig^it) + Q{S{t) - g^it)). (6) 

Whereas the first term on the right hand side describes the active relaxation of the tagged 
subsystem, starting in state z, the second term characterizes a passive exchange process from 
some state j (j 7^ i) to state i due to the coupling to the adjacent subsystem. Since we aim 
to keep the dependence on Q only up to linear order, we can substitute S{t) by S^(t). Then 
integration yields 

-l = -(r, + Q)G, + Qr°. (7) 

Thus, one can write 



Multiplication with pi and summation over all states yields 

-" = -°-Q^((:^a/r^)-i)- (9) 

Note that (l/F^) is identical to J dttS^{t). It is exactly the first term in the brackets which 



m 



27| has been denoted as 1/Pm and is a measure for the non-exponentiality of the total 



relaxation, i.e. of the dynamic heteterogeneities. Thus, we can rewrite 

More generally, in lowest order of q the tagged system of Nmm particles will experience a 
coupling to the surrounding particles which is proportional to the number of particles, i.e. 

— Nmin- Furthermore, after identification of r° with Ta{Nmin) one may finally write for 
N > Nfnin to lowest order in q 

r.(iV) ^ J N \ T^{N^,^) / I \ ^^^^ 



TaiNmin) V^mm / (t) \I3m J 

Thus, after determination of /3m and TaiJ^min)! this relation can be used to extract 
the dimensionless coupling constant q. 

Note that this estimation breaks down if (i) the second term on the right hand side 
approaches one (then higher-order terms of q matter) or (ii) starts to be much larger than 
Nmin- Then the coupling is no longer proportional to because the variation of distances 
between different subsystems starts to matter. 



16 



